Anderson Lattice Description of Photoassociation in an Optical Lattice 
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We consider atomic mixtures of bosons and two-component fermions in an optical lattice poten- 
tial. We show that if the bosons are in a Mott-insulator state with precisely one atom per lattice, the 
photoassociation of bosonic and fermionic atoms into heteronuclear fermionic molecules is described 
by the Anderson Lattice Model. We determine the ground state properties of an inhomogeneous 
version of that model in the strong atom-molecule coupling regime, including an additional har- 
monic trap potential. Various spatial structures arise from the interplay between the atom-molecule 
correlations and the confining potential. Perturbation theory with respect to the tunneling coupling 
between fermionic atoms shows that anti-ferromagnetic correlations develop around a spin-singlet 
core of fermionic atoms and molecules. 

PACS numbers: 03.75.Ss, 05.30.Fk, 32.80Pj, 67.60.-g 
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Ultracold atoms and molecules trapped in optical lat- 
tices provide an exciting new tool for the study of 
strongly correlated many-body systems 0] . The exquisite 
degree of control of the system parameters permits the 
detailed study of a variety of exotic states of matter, and 
as a result these systems are contributing to the estab- 
lishment of significant new bridges and interplay between 
AMO science and condensed matter physics. While much 
work along these lines has concentrated so far on ultra- 
cold atoms 01 , the coherent formation of bosonic or 
fermionic molecules 0: IM S via either Feshbach reso- 
nances 8j or two-photon Raman photoassociation 9] of- 
fers an additional path to the study of strongly correlated 
atoms and molecules ^3 Very recently, collective 
coherent phenomena between an atomic and a molecular 
gas in an optical lattice have been observed experimen- 
tally 0. 

In this letter we analyze the ground state of a mixture 
of atomic bosons and two-component fermions coupled 
to heteronuclear fermionic molecules [l^lll^ by photoas- 
sociation or Feshbach resonance in an optical lattice. We 
show that this system can be mapped onto the Ander- 
son Lattice Model (ALM), a model that has previously 
found important applications in the description of heavy 
electrons and intermediate valence systems in condensed 
matter physics. In particular, this model is known to ex- 
hibit a great variety of possible behaviors, such as e.g. 
the Kondo effect and magnetic ordering |l5| . 

In the context of AMO experiments, the modifications 
of the ground-state properties due to the presence of a 
trapping potential are of particular interest. For exam- 
ple, in the bosonic (fermionic) Hubbard model, that po- 
tential is known to result in the coexistence of Mott- 
insulator and superfluid (metal) phases [lil Il7| . In the 
strong atom-molecule coupling regime under considera- 
tion here, the ALM exhibits two types of magnetically 
correlated states of fermionic spins on the lattice, an 
on-site spin-singlet (paramagnetic) state of the atoms 
and molecules, and an anti-ferromagnetic (AF) corre- 



lated state among lattice fermionic spins. The inhomo- 
geneity of the confining potential gives rise to a spatial 
structure of coexisting paramagnetic and AF domains. 

We consider at zero temperature a mixture of atomic 
bosons of mass and atomic fermions of mass Ala 
and spin a trapped in an optical lattice poten- 

tial The fermionic and bosonic atoms can be coher- 
ently combined into heteronuclear fermionic molecules of 
mass Aim = A'h -t- Ala via two-photon Raman photoasso- 
ciation. The lattice lasers are tuned so that the fermionic 
atoms experience a weaker on-site lattice trapping fre- 
quency, r^Q, than the bosonic atoms and the fermionic 
molecules, i7fc,r2,„ ^ ^a- 

The Hubbard-type Hamiltonian describing this system 
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Here aia , rhiu , and hi are the annihilation operators of the 
fermionic atoms and molecules and of the bosonic atoms 
at the z-th site, respectively. The corresponding number 
operators n°;„ = a\^aia, h™. = ml^rhia, and n- = blbi 
have eigenvalues n°^, n^, and n^. 

The on-site single-particle energies e^, where a = 
a,b,m, consist of the sum of the contributions e" ~ Qa 
from the lattice and V^"^ from the trap potential at site i, 
the two-photon detuning between the Raman-lasers and 
the difference in internal energies between the atoms and 
molecules being included in ly. In case the external poten- 
tial is created by optical means, the spin structure of the 
fermions is not resolved, resulting in spin-independent 
single particle energies. This is the case that we con- 
sider here. The terms proportional to Ua, Um, and Ub 
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describe on-site two-body interactions, inter-species colli- 
sions being ignored for simplicity. Finally, the interaction 
Hamiltonians Hn describes the conversion of fermionic 
and bosonic atoms into fermionic molecules, and tun- 
neling between nearest neighbor sites denoted by {ij) is 
described by the parameters ta, t„i, and if,. 

We concentrate on the case Uh,m ^ g ^ h^m and 
g ^ ta ^ Ua, a regime that is achieved for e'''™ ^ e" 
and Mb 3> Ma- One can then ignore intersite tunnel- 
ing of the fermionic molecules and bosonic atoms, whose 
rate is given by the exponential of the ratio between the 
lattice constant A/2 and the width of the localized state 
la = {^'^ /MaSa)^^'^ 11]. We assumc the opposite con- 
dition for the fermionic atoms and neglect instead their 
on-site two-body collisions 22]. Typical parameters for 
^^Rb atoms are et/h ~ 10^ s^^ and Ub/h ~ 10"* s~\ and 
the photoassociation coupling constant is estimated to 
be of the order oig/h~ 10^ s^^ 0. Since heteronuclear 
diatomic molecules are polar, we also need an estimate 
of their electric dipole-dipole interaction. The dipolc mo- 
ment of the triplet state of the KRb dimer is 0.02 
a.u. 12^, (1 a.u. is 8.478 x lO'^" Cm), in which case 
the dipole-dipole interaction energy is less than 10 s~^ 
for A — 1064 nm and can safely be ignored. We further 
neglect all loss processes. 

In the absence of photoassociation, the ground state 
of the atomic fermions consists of free fermions in the 
confining trap V"" while for bosons it is a Mott-insulator 
core with exactly one boson per site within a set {f\ of 
sites from the center of trap up to edges determined by 
the number of bosons Nb, and zero occupation outside. 
We therefore have 

Y,fC + n\ = liovze{i}; J2M.^Nf, (4) 

where n£. fi^^ + n™. The former constraint is due 
to prohibition of double occupations and neglect of tun- 
nelings for bosonic atoms and fermionic molecules |23| . 
These restrictions result in the reduced Hamiltonian 

Hr = -ta J2 + H-C.) + J2 (5) 

+ I^l'^ + er - +gY,{rntka^, + H.c). 

ia ia 

In the homogeneous case Vf = 0, this Hamiltonian is 
equivalent to the so-called slave boson model, which un- 
der the constraints (0)) reduces to the periodic ALM with 
infinite on-site repulsion. We note, however, that in our 
model the bosonic operator b corresponds to a real phys- 
ical particle, in contrast to the pseudo-boson introduced 
in the study of the Kondo problem . In the following 
we investigate the ground state properties of this system 
in the presence of a harmonic trap and for weak tun- 
neling of the fermionic atoms. Since the atom-molecule 
(AM) correlation energy is order of g as shown later, this 



is justified in the regime g ^ ta- We proceed by first ap- 
plying a counting argument for the case ta = Q under the 
constraints and then account for tunneling coupling 
in perturbation theory. 

For ta = 0, the ground state is a product of lowest 
energy states with total fermion number = n°' + 
_ J2^n°;^ + J2a''^ia Bach. lattice site, l$o) = 
Hi l?^''^; lowest energy) i. For lattice sites belonging to the 
set {£}, we can have Uf = 0, . . . , 3 fermions, and the 
corresponding normalized lowest energy states are 

\nf =0)l = \0a,0rn)l, (6) 
|r/ = l;a)i = ai\aa,Om)l + Pl\Oa, CTm)l, (7) 
Jr/ ^ 2)i ^ Ui\2a,0„,)l+Vi (I ta,ir«>( - | ia,Tr«>;X?) 

\nf ^3;a)i = \2a,am)i, (9) 
with corresponding energies 
Ei{0) - 0, 

Ei{l) = et + 5i-^Sf+g^, 

Ei{2) = 2e« + 6i- ^6f + 2g\ 

Ei{3) = ie1 + 25i, (10) 

where we have used the compact notations 
\n(Ta, n'(7'^)i = \n(7a)i ® \n'a'^)i ® jn'' = 1 - n')i and 2Si = 
ly+ef-e'l-ef. Here a/ = -g{[{Sf+gY^^-^i]^+g^y^^^ 
and ui = ~V2g{[{Sf + 2g2)i/2 _ Si]^ + 2g^y^/^. The 
fermionic atom and molecule populations at a single site 
are {nf) = af and (n™) = Pf respectively for the — 1 
manifold, and (nf ) = 1 + < 2 and (nj") = 2vf < 1 
for = 2. We observe that the average number of 
fermionic atoms gradually decreases with decreasing 
(5;, while the number of molecules increases. From 
Eqs. 10 and (O, we also note that the atom- molecule 
conversion leads to correlated states in the n-^ = 1 and 
= 2 manifolds, with energies lowered by an amount 
of the order of g. The state jn-^ = 2) corresponds to 
a spin-singlet state, while spin-triplet states such as 
ferromagnetic configuration have no coherent scattering 
and then are higher energy states. 

To count the number of fermions and find the ground 
state configuration, we introduce at each lattice site I 
the local chemical potential r]i ( M) = Ei (M) -Ei{M-l) 
{M — 1, 2, 3). It corresponds to the energy needed to add 
one fermion to M — 1 fermions already at that site. At 
the sites i ^ {£}, the local chemical potential is equivalent 
to the atomic single particle energy, -qj — rjf = e°. We 
then proceed by first determining the Fermi energy ep 
for a fixed number of fermions, the fermionic occupation 
number at a given site being then determined by the 
condition rjiin^) < ep l£ '7!("-'^ + l)) where we set rii{0) — > 
— oo and t?/(4) +oo for convenience. 

In order to gain insight into the ground-state properties 
of the system, we consider the specific case of a one- 
dimensional, harmonic trap potential V°' for each species 
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FIG. 1: (Color online) ryi(M) for S/g = and Ni = 64, 
and for (a) V/g = 1, (b) V/g = 2.5; Ground state density 
distribution of (c) fermionic atoms, (d) molecules for 5/g = 0, 
V/g = 1 (solid line), S/g = 0, V/g = 2.5 (dashed line), and 
5/g = -1.5, V/g = 1 (dot-dashed line) for Nf = 2Nb = 128. 

a = a, 6, m [13, = F"{2/(iVfc- l)}2{i- (AT, + 1)/2}2, 
with the trap center located half-way between two lattice 
sites. We assume that the strength of molecule trapping 
is the sum of those of the atomic fermions and bosons, 
i.e., = V"" + , resulting in the uniform detuning 
25i = 25 = V + e™ — e" — e''. We also set the offset energy 
= without loss of generality. 

Figures ^,b show the local chemical potentials ?7;(M) 
at each lattice site for (5 = and for (a) — 1 

and (b) V"" / g — 2.5. In this example, the system would 
contain 64 bosonic atoms occupying lattice sites within 
1 < Z < 64 in the absence of photoassociation, g = 0. 
We remark that in that limit all rji{M) become degener- 
ate and equal to e° at (5 = 0, clearly indicating that the 
gaps among the three chemical potential "bands" arise 
from the AM coupling. The corresponding Fermi ener- 
gies for Nf = 2iVf, = 128 are the solid lines. In case (a), 
the ground state corresponds to a configuration where all 
sites within {C\ are in the spin-singlet state \nf — 2), all 
other sites being in the = 0) state. In contrast, in 
case (b) the trap energy near the edges of the atomic den- 
sity overcomes the AM correlation energy at the center 
of trap. As a result, the fermions are in the uncorrelated 
state = 3) near the center of the trap and in the 
correlated state = 1) at the edges. 

The ground-state density distributions of fermionic 
atoms (n") and molecules (n™) are shown in Fig.^J; and 
respectively. For V°-/g — 1 (solid fine), spatially ho- 
mogeneous densities of fermions are realized up to the 
edge of the sites {1} and form the spin-singlet state of 
Eq. 0. This corresponds to the Kondo insulator 
state in homogeneous systems. Because of the correlated 
gap energy this state is stable against small amounts of 



atomic tunneling ta- For V"" /g = 2.5 (dashed line), atoms 
and molecules form a spatial shell structure reminiscent 
of the Mott-insulator shells of bosons or fermions in com- 
bined optical lattices and trapping potentials 

As 

5/g becomes large and negative, the molecule popula- 
tion at each lattice site approaches unity, but the inho- 
mogeneous character of the population persists even for 
V"' / g = 1 at 5/g = —1.5 as shown by the dot-dashed 
lines in Fig. ^,A. 

We now discuss the effects of weak intersite tunnel- 
ing of the fermionic atoms. The role of the tunneling 
is two folds. First, the enhancement of fluctuations of 
fermion number on the border between different shells. 
Second, the anti-ferromagnetic (AF) correlation on the 
\nf = I] a) domain. In the strong AM couphng regime, 
g ^ ta and |^| < g, these can be investigated by (quasi- 
) degenerate perturbation theory up to second order in ta- 
The ground state is given by 

l*9>- E Cg{n{,a,}\Mn{,a,}). 

The set of {n{,ai} contains those fermion number con- 
figurations, {n{}, which can be realized through multiple 
tunneling with an only energy cost {V"' — T^^^K^C g) for 
each process from the unperturbed ground state in the 
absence of ta- A possible spin configuration {ui} can 
couple to a different one via virtual excitations of energy 
lager by an amount of order g. 

To illustrate the role of the number fluctuations we 
consider the case Nt = 16, V/g — 1.5, and 5 — with 
Nf = 22 such that the unperturbed ground state in the 
system of = with \n-^ = 1) domains surrounding the 
spin-singlet core oi \nf = 2) and no \n-f — 3) compo- 
nent. Figure |2] shows the density and number fluctua- 
tions of fermionic atoms, defined by (erf )^ = {{{n°;Y) — 
{fif)^) / {fif) , as a function of ta/g for the direct diago- 
nalization in the first-order perturbation theory. With 
increasing ta the hybridization of the different n-^ states 
is higher on the border between = 1 and = 2 
shells, while the spin-singlet core and \n^ = 1) domains 
remain. Consistently with this result, the number fiuc- 
tuations of atoms develop strongly on around the border 
region and weakly on the core and surrounding regions. 
The reason of the suppression of fluctuations on the edge 
of the sites {1} is that atomic tunneling from the edge 
sites into outside costs an energy ~ g, so that such a 
transition is neglected. We note that the result of Fig. [3 
does not depend on qualitatively the exeat value of the 
detuning parmeter and number fluctuations of molecules 
are vanishingly small for any \5\ < g as expected. 

The second-order corrections of atomic tunneling be- 
tween nearest neighbor sites that are both within the 
\n^ — 1) manifold leads to an AF correlation among 
fermionic spins composed of mixtures of atoms and 
molecules. In contrast, if at least one site is in a n-'' 7^ 1 



4 




FIG. 2: (Color online) (a) (n") and (b) erf as a function of 
ta/g for Nb = 16, Nf = 22 and for = 0, = 1.5. 



manifold, those tunneling processes are spin-independent 
and result only in an energy shift. The origin of AF corre- 
lation is similar to that observed in Hubbard model in the 
strong coupling limit '2l'|, while in our problem at hand 
the coherence character between atoms and molecules 
plays a key role. 

Introducing the coherent spin operators at the i-th site 



J 



(l/2)(4-nf^). 



/ ^ J 



the spin-dependent effective Hamiltonian is 



f i 



(11) 



As a consequence of AM coherence, the AF coupling coef- 
ficient Jy (> 0) depends on the detuning parameter and, 
when I V^" — ^ g, Jij is reduced to a site- independent 
one J = t^/g • js(<5) shown in Fig. 0| The last term of 
the above Hamiltonian is owing to possible tunneling be- 
tween a site occupied by an atom and a neighboring one 
occupied by a molecule even for parallel spins. 

Since Hamiltonian Ijlll) is spin 1/2 Heisenberg model, 
we can make out it's magnetic properties on a — 1) 
domain. If we consider a one-dimensional example like 
Fig. 12 Hamiltonian Hll|l affects left and right sides sep- 
arated by spin-singlet core independently. The ground 
state corresponding to each domain is characterized by 
^ 0- The asymptotic form of the spin den- 
sity correlation function has a power law behavior 

(sj+fc jsj^z) ^ ("l)''/^ except for finite size corrections. 

The present mechanism of AF correlations is unique in 
that atomic tunneling and strong AM couplings act coop- 
eratively, resulting in coehrent AF characteristics among 
intersite fermionic spins of both atoms and molecules. 
Moreover, the AF correlation develops outside the para- 
magnetic spin-singlet core, so that the spatial param- 
agnetic and AF correlated domains coexist. For find- 
ing magnetic properties on the border region or when 
J ^ ta ^ g, need to consider the interplay of number 
fluctuations and AF correlations. This will be discussed 
in a future publication. 

In summary, we have studied atomic boson-fermion 
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FIG. 3: Scaled AF coupling strength js as a function of 5. 



mixtures with atom-molecule conversion in the presence 
of an optical lattice potential plus a harmonic oscilla- 
tor potential. We have shown that by controlling lat- 
tice depths for fermions and bosons independently, an 
inhomogeneous version of the ALM can be realized. The 
new feature of inhomogeneity gives rise to various spa- 
tial structures of correlated atom-molecule ground states. 
We have analyzed the tunneling coupling effect by de- 
generate perturbation theory and have shown that this 
coupling and inhomogeneity give rise to coexistence of 
spatial domains of paramagnetic and AF correlations. 
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